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Abstract. We show how the combined use of the generating function method 
and of the theory of multivariable Hermite polynomials is naturally suited to 
evaluate integrals of gaussian functions and of multiple products of Hermite 
polynomials. 



The generating function method is often exploited to derive the analytic form of 
integrals of the type [1] 

/oo 
dxH n (ax + b,y)e~ cx2+ax 
-oo 

where 

I"/ 2 ! T n-2fc„.fc 

is the two variable Hermite polynomial with generating function [2] 

1 — ■ - -r.t+yt 2 



(3) Y.7A H ^y) = eXi 

L — ' n! 

n=0 

By taking into account this identity, from eq. ([T]) we obtain 

00 j-n poo 

(4) V L r / n =e i,t + y * 2 / dxe- cx2+ ^ at+a ^ , 

nl J -oo 



n=0 



and, thus, we have reduced our problem to the evaluation of a trivial gaussian 
integral, which yields 



n=0 V k 



a 2 /, a a 



(5) yy-in=^ex P {f- c +[b+f^) t+[y+f- c ) t 2 



The use of the generating function allows us to write the integral ([T]) as 
follows 

(6) 7 « = ^^jM^'^ 
Moreover, since 

(7) H n (x,y)=exp(yd 2 x ) x n , 
we can write I n in the following operational form 

(8) 
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Let us now consider the integral 

/oo 
dxx m H n (ax + b,y)e- cx2+ax , 
-oo 

that occurs in some problems involving quantum harmonic oscillator, or, in classical 
optics, in the evaluation of overlapping of Gauss-Hermite beams. It can be cast in 
the form 

(10) m /„ = d™ I n , 
from which, taking into account the identities 

(11) d*H n {x, y ) = T ^-rr i H n - k {x,v) dle^ 2 = H k (2 0x,0) e' 3 ** , 

(n — k)\ 



one 



finds 



Jtt I a 2 \ ( a 1 a a a 2 a \ 

(12 m I n = —= exp — H m , n [—,—;b+—,y + — — , 

y c \ 4 c/ \ 2 c 4 c 2c 4c 2 c/ 

where we have introduced the two-index Hermite polynomials [U [3] defined as 

follows 

(13) 

min (m,n) 

H m ,n{x,y,w,z | r) = — _ — - — - r k H m . k (x, y) H n . k (w, z) . 

The relevant generating function writes 

(14) } — - — - H m n (x, y : w, z I t) = exp (:cu + ?/it 2 + i;u; + ;z?j 2 +Titi;) , 
£ — ' m! n! 

m,n— 

and, thus, using the same procedure as before, we find for integrals of the type 

/oo 
dx H m (a x + b, y) H n (c x + d, z) e ~ s x2+a x 
-oo 

the explicit form 



7T / a 2 



(16) im,n = —j exp ^— J ff m ,„ (x, y; w, z I r) , 
where we put 

a _ a 2 _ c _ c 2 a c 

x = b H a , v = V -\ , iu = a H a , z = z H . r = — . 

2/ ' y 4/' 2/ ' 4/' 2f 

By taking into account that 

(17) H m>n (x,y;w,z | r) = exp (y d% + w d 2 + t d 2 z ) x m y n 
we obtain the generalization of eq. (JSJ as 

(18) I m , n = ^ exp (^j exp (y d 2 + w d 2 + r d 2 s ) x m y n . 

The operatorial method can be applied from the very beginning and the integral 
(TTBl) can be written as 



(19) I m , n = exp (yd 2 + z d 2 ) l m ^ n , 
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oo 

I m , B = / dx(ax + b) m (cx + d) n e~ fx2+ax 

— oo 

/tt ( a 2 \ „ / a 2 c 2 , 



where 



(20) = ^=exp^— j ff m , n ^,— ;tfl,— |r 

This approach may simplify the computation of integrals involving products of more 
than two Hermite polynomials. 

It can now be shown that for the integral 

/oo 
dxx p H m (ax + b,y)H n {cx + d,z)e- fx2+9X , 
-OO 



an identity analogous to the (|T0|) holds 

(22) plm.n — $a. ^Tn,Ti • 

The use of the identities 

Tfl\ 

d x H m n (x,y,w,z | t) = — - H m - k ,n(x, y;w,z | r) 

(to — ky. 

n\ 

(23) d w H mn (x,y;w,z | r) = — - H m ^ k {x, y; w, z | r) , 

[n — ky. 



allows us to obtain for the integral (f21[) the following expression 



where 

= d*H m , n {x,y;w,z\ r) 



, k k 

c 



( 25 ) = { f] ) g Q S (TO-/)!(n-fc + 0! Hm ~ j ' n - k+l ( "' R ®" " ' T) 



& to.' n) 



As a final example, we consider the integral 

r°° (ax + b) n 
(26) I n = j_Jx [ {l+ + c J r (n<2u,ueR). 

The combined use of the generating function and Laplace transform methods 

yieldfO 



oo ,„ bt 



(27) £ = / cUe-'*"- 1 / dxe- 

which can be finally explicitly worked out as 



sex +a x t 



(28) I n = -£— H%> [b,— 

VcT(v) \ 4c 



In performing the sum we do not take into account of the condition n < 2 u, but the consequent 
result in eq. Q28|l must be considered valid only under this hyphotcsis. 
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where 

[n/2] tv u 1 

(29) ^)-^£V.-^»«-'V. 

These polynomials do not belong to the Appell family and the relevant generating 
function can be written adopting an umbral notation as follows 

(30) Y,- H n ) ^y)= eXt+1)t y k = T(v-k-l/2)y k . 

n=0 U - 

With this assumption we also get 

(31) Hi»\x,y)=exp(yd 2 x ) x n , 

and then most of the properties of this family of polynomials can be derived straight- 
forwardly. 
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